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0. INTRODUCTION 

Among the natural variational integrals in geometry are the inevitable integrals on space 

curves 7(5); these include length £(7) = / ds, total torsion T(j) = f rds, total squared 

curvature £(7) = \ J K 2 ds, the integral ^(7) = h J K 2 Tds — i j j s x 7 ss -7 SS s ds, and an 

infinite sequence of integrals following these four. A closely related sequence, for curves 

T(a) on the sphere, begins with total geodesic curvature A(T) = f nxdsx, and energy 
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£{T) = J \\T a \ 2 da. (By the Gauss-Bonnet theorem we can view A(T) as the spherical 
area bounded by T, when T is closed). Though the whimsical term inevitable will be 
explained here by some simple geometric considerations, these sequences first arose via 
some physical models (which the geometric point of view has helped illuminate). 

In fact, the length integral £ may be regarded as the Hamiltonian for a model of 
thin vortex tubes in three dimensional hydrodynamics, where the evolution of the tube's 
centerline 7(s, t) is governed by the Localized Induction Equation (LIE): 7* = 7 S x 7 SS = 
kB (see [M-W]). Then the remaining integrals in the first list are constants of motion 
for the evolving curves. Similarly, the energy £ is the Hamiltonian for the (continuous) 
Heisenberg Chain (HC) (see [F-T]), T t = T x T aa , while A and the integrals following £ 
are constants of motion for the curves T(cr, t) evolving on the sphere. 

The simple geometric relationship between these two soliton equations is this : if 7 
satisfies (LIE), then the tangent indicatrix of 7, T = 7 S , satisfies (HC). By the way, 
it pays not to check this fact too hastily. After differentiating the first equation with 
respect to s, one wants to replace "fts by 7 s t = T t to get the left hand side of (HC). To 
justify this one has to observe that W = kB happens to belong to the special class of 
locally arclength preserving (LAP) variation vectorfields. W is LAP if the speed along 7 is 
preserved under the evolution 7^ = W - in particular, an arclength parameter s remains 
an arclength parameter - which is equivalent to the condition that W s is orthogonal to 
T (obviously satisfied by W = kB = T x T s ). 
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The two lists of constants are also related by the tangent indicatrix construction; if 
T = 7 S , then T(i) = A{T),K{i) = £(T), etc. (Note that the second list is 'missing' 
one integral, namely C.) The two lists are essentially 'equivalent', but the first is of 
perhaps greater interest in the study of geometric variational problems. For example, 
one of the oldest problems in the calculus of variations is that of the Bernoulli elastic 
rod, to 'minimize' /C with a constraint on C. More generally, the extrema for sums 
a±C + a^T + a^/C, ctj constants, are precisely the (centerlines of) Kirchhoff elastic rods 
in equilibrium. (Actually, this simple characterization seems to be relatively unknown; 
the twist energy, which the Kirchhoff model adds to the Bernoulli model is given by a 
quadratic integral on framed curves [L-S]). To our knowledge, no one has classified the 
extrema for sums a±C + C12T + CI3/C + a^T . 

The first of these variational problems was introduced into the study of (LIE) in 
1971 by Hasimoto, who observed that if 7(3, 0) represents a (Bernoulli) elastic rod in 
equilibrium, then 7 moves through space under (LIE) without changing shape [Has 1]. 
More generally (not part of Hasimoto's original observation) such one-soliton solutions 
to (LIE) are precisely the equilibria for the Kirchhoff model. 

With this beginning, Hasimoto went on to relate the filament equation itself to yet 
another physical model - already known to exhibit soliton behavior - the Non-linear 
Schrddinger Equation (NLS) : —i^t = (i> ss + ||^| 2 ^)- Hasimoto showed that if a space 
curve 7 evolves according to (LIE), then the curvature and torsion of 7, k(s, t) and r(s, t), 
can be combined into a complex curvature function = «'/•*, which satisfies 
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(NLS) [Has 2]. The proof of this fact amounts to computing variational formulas for 
k and r or equivalently, for the complex curvature ip. It turns out that these formulas 
have a special structure, which is related to the infinite list of constants of motion for 
(LIE) , and which also clarifies the precise nature of the equivalence of the three models. 
The study of these formulas and their relatives, their geometric interpretation, and their 
geometric applications form the central theme of this paper. 

To describe this special structure, let ^(7) = ip denote the Hasimoto transformation 
taking space curves to complex curvature functions. Then the desired variational formula 
for tp is contained in the following general formula for the differential of H at 7 in the 
direction of W , a LAP vectorfield: 

(1) dH(W) = -zn 2 (w). 

The most important feature of this formula (whose proof is included below) is the ap- 
pearance of the (squared) recursion operator 71, a linear integro-differential operator on 
vectorfields. The operator Z can be thought of as a simple isomorphism from LAP 
vectorfields to complex functions of s; formulas for 1Z and Z will be given below. We 
write = to denote equivalence modulo addition of terms of the form icvp = ic7i,(j),c 
a real constant. Note we have defined H itself only up to multiplicative factors of the 
form e* 9 , 6 a 'constant of integration' - this ambiguity is not reflected in k or r. For a 
varying curve, differentiation of 9 = 9(t) produces an extra term icip in the variation of 
ip. Setting W = kB in (1) easily yields ip t = dH{W) = i(ip ss + ±(|^| 2 + c(t)))ij>), which 
is the precise statement of Hasimoto's result. 
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1. THE RECURSION OPERATOR AND NATURAL FRAMES 

To discuss formulas for 1Z and Z we use natural frames along 7, in place of the standard 

Frenet frame T = {T, N, B}. Recall the Frenet equations may be written in the form 

T s = QxT,N s = QxN,B s = QxB, where O is the Darboux vector (s- angular velocity 

of J 7 ) given by uj = tT + kB. Similarly, the natural frames are those orthonormal 

frames Af = {T, U, V} along 7 having s-angular velocity u s = —vll + uV ', for some 

'curvature functions' u = u(s) and v = v(s); so Af satisfies the 'Frenet Equations' T s = 

uj s xT = uU + vV, U s = -uT, V s = -vT . Note that u s = T x T s = kB, the vectorfield 

defining (LIE)! In other words, among all (adapted orthonormal) frames N = {T, U, V}, 

the natural frames are characterized by: M has zero tangential component of angular 

velocity . One might prefer to call M natural if it has constant T-component of angular 

velocity - think of a bead sliding and turning without friction along a wire 7 - but we 

will not use this more general class here. 

An alternate way to describe natural frames is to say that the vectors {U, V} normal to 
7 are obtained by parallel translation in the sphere along the tangent indicatrix T. Note 
that while T is uniquely determined, Af is determined only after an 'initial frame' has 
been chosen at some point along T. An advantage of M (which plays no role presently) 
is that Af can be defined even where 7 has vanishing curvature. 

We give formulas for Z(W) and 1Z(W), where W = aT + bU + cV, and we also provide 
a simpler expression for the complex curvature function ip\ 
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ip = W(7) = U + iv, 
Z{W) = b + ic, 

(2) 

n(w) = -v(Txw a ), 

V(W) = { J bu + cv ds)T + bU + cV 

Here we have also introduced a parameterization operator V, which leaves the nor- 
mal part of an arbitrary vectorfield W alone but turns W into a LAP vectorfield (since 
j^V(W) has no T component). It follows that the recursion operator 1Z preserves the 
class of LAP vectorfields. Note that we are again being casual about constants of inte- 
gration; V(W) has only been defined up to addition of terms of the form constant ■ T, 
and by the same token, Z is not quite an isomorphism of LAP vectorfields. 

Until now, we haven't said much about the meaning of the recursion operator 1Z - 
we have only indicated that it allows for a very compact expression for the variation of 
the complex curvature. To get some experience with 7Z, let's apply it to the humblest 
of all LAP vectorfields: 1Z(T) = -V(T x T 8 ) = -V(kB) = -kB, an old friend! In this 
computation we simply ignored V since kB is already LAP. 

Since this worked out so nicely, let's apply 1Z again (keeping in mind that tangen- 
tial terms inside of V{ ) can be discarded): TZ 2 {T) = K(-kB) = V(T x (T x T s ) s ) = 
V(-T ss ) = -V(u s U + v s V) = -(J u s u + v s vds)T-u s U-v s V = -\{u 2 + v 2 )T -u s U - 
v s V. Note our good fortune in being able to compute the antiderivative explicitly! As 
it turns out, this will continue to happen as we successively apply 71, and the LAP 
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vectorfields X n = 1Z n (—T), n=l,2,3,..., are none other than the Hamiltonian vector- 
fields of the constants of motion for the Localized Induction Equation - hence the names 
"recursion operator" and "inevitable integrals" . 

2. THE VARIATIONAL FORMULAS 

We have explained the name, but only part of the significance of 1Z. Recall that the 
formula (1) relates variations in the 'first' model (LIE) to variations in the 'third' model 
(NLS). Since the square of 1Z occurs in this formula, it seems reasonable that a first power 
of 1Z ought to relate the 'second' model (HC) - this is clearly the intermediate model - to 
(LIE) and (NLS). To see how this is true, it is useful to modify the (HC) model slightly, 
replacing the curve T by a natural frame Af = {T, U, V}. Note that in doing so, we have 
only enlarged our 'state space' by one dimension, corresponding to our usual ambiguity 
(say, the choice of 'initial frame' {Uo, Vb}, or choice of 'phase factor' e 10 for ip). 

Now suppose Af(s, t) evolves in time according to T t = u l x T, U t = uj 1 x U, V t = 
u l x V, for some t-angular velocity u) 1 = u/(s, t). For Af to represent an evolving natural 
frame, it must satisfy = U s ■ V, for all s, t. It follows that for u l to be naturality- 
p 'reserving (NP) it must satisfy: = (U s • V) t = U ts • V + U s ■ V t = (u* x U) s ■ V + 
U s ■ (uj 1 x V) = (o/) s x U ■ V = (uj t ) s ■ T, a familar condition! In other words, allowing 
the same vectorfield (aT + bU + cV) to play two different roles: W = aT + bU + cV is 
LAP if and only if u f = aT + bU + cV is NP. In particular, this means that if u/ is NP, 
so is n n (uj t ),n= 1,2,3, .. .. 
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Without further ado, we state the promised relationship between the variations of 
curves and variations of the corresponding natural frames: If W is LAP, and j(s, t) 
evolves according to j t = W, then the t-angular velocity of a natural frame Af = {T, U, V} 
along 7 is given by 

(3) u* = -nw 

Proof: Since W = aT + bU + cV is assumed to be LAP, we have 7^ = 7^, so tu t x T = 
T t = 7 ts = (aT + bU + cV) s = (au + b s )U + (av + c s )V . It follows that u* = rfT - 
(av + c s )U + (au + b s )V, for some function n = r)(s,t). To determine 77, we compare 
expressions for U st and U ts : U st = (—uT) t = —u t T — u(au + b s )U — u(av + c s )V, 
while U ts = (u l x U) a = ((r]T + (au + b s )V) x U) s = (-(au + b s )T + V V) S = -(b ss + 
a s u + au s + rjv)T — u(au + b s )U + (r) s — auv — vb s )V ', so equating "^-coefficients gives 
i] s — b s v = —uc s . Therefore u l = (f vb s — uc s ds)T — (av + c s )U + (au + b s )V. Finally, 
-K(W) = V(T x (aT + bU + cV) s ) = V(T x ((au + b s )U + (av + c s )V) = V(-(av + 
c s )U + (au + b s )V) = (f -(av + c s )u + (au + b s )vds)T - (av + c s )U + (au + b s )V = w*. 

It remains to see how 1Z takes us from the second model to the third model. Given 
a complex function ip = ip(s,t), we can always write the t-derivative ip t in the form 
ipt = — -2(0*)) for some vector field 0* = aT + bU + cV . Of course, 0* is not unique since 
Z ignores tangential terms. But if 0* is also required to be LAP, then it willbe unique up 
to the usual terms, constant ■ T. Now the question is, if a natural frame M = {T, U, V} 
evolves with t-angular velocity a/, and if = u + iv = Z(ull + vV) = Z(T S ) evolves 
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accordingly, what is the vectorfield 0*? The required computation is simplified by the 
abuse of notation Z = U + iV, Z(W) = Z -W (complex inner product): ip t = (Z ■ T s ) t = 
Z t -T s + Z-T ts = iv* xZ-T s +Z-(u f xT) s = Z ■ (-T x = Z(ft(u;*)). In 

other words, we have the following formula relating the variations in natural frames and 
complex wave functions: 

(4) = -m u *). 

Finally, note that in case J\f(s,t) comes from a curve j(s,t), and ^y t = W is LAP, we 
have ip t = ZiTZ^)) = Z(K(-K(W)) = -ZTZ 2 (W), which proves (1). 

In order to summarize the various roles of the recursion operator 1Z we have included 
a diagram below. Note that 1Z connects levels corresponding to the three models (LIE), 
(HC), and (NLS) (for the actual equations (LIE), (HC), and (NLS), set n = 1 ) and K 
also connects t and s directions of motion. Some of the arrows have not been explained, 
but are easily interpreted and checked. 



Is > It 



-n 



s -K n t 

u s > ur 



-n 



-n 



We have not discussed here the Hamiltonian nature of the equations (LIE), (HC), 
or (NLS), or the actual definitions of the spaces on which these Hamiltonian flows are 
defined. For background on such points see [L-P 2], where (1) is first proved and used to 
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establish the precise nature of the Hasimoto transformation as a Poisson map. Arguments 
similar to those used in [L-P 2] can be combined with formulas (3) and (4) to give 
similar characterizations of the intermediate maps discussed here (thus interpreting Tt 
as a composite of Poisson maps). However, our intent here has been to emphasize the 
ubiquity of the recursion operator 1Z and its relation to geometric varational problems 
and flows of geometric origin. 
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